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Some aspects of electronic topological transition in 2D system 
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oo ■ 

(N| ; M.Kiselev 1 ' 2 ' 



F. Bonis 1 , F.Onufrieva 1 and P.Pfeuty 1 



S3 
O 
O 

Oh' 



o 

> 

O 

00 

On 

-i— > 
c3 



1 Laboratoire Leon Brillouin, CE-Saclay 91191, Gif-sur-Yvette, France 

2 Russian Research Center "Kurchatov Institute", 123 182 Moscow, Russia 



I 

C 

O 

o 



X 



February 1, 2008 

Abstract. We study the ordered "excitonic" states which develop around the quantum critical point (QCP) 
associated with the electronic topological transition (ETT) in a 2D electron system on a square lattice. We 
consider the case of hopping beyond nearest neighbors when ETT has an unusual character. We show that 
the amplitude of the order parameter (OP) and of the gap in the electron spectrum increase with increasing 
the distance from the QCP, 8 C — S, where 5 — 1 — n and n is an electron concentration. Such a behavior is 
different from the ordinary case when OP and the gap decrease when going away from the point which is a 
motor for instability. The gap opens at "hot spots" and extends untill the saddle points (SP) whatever is 
the doping concentration. The spectrum gets a characteristic flat shape as a result of hybrydization effect 
in the vicinity of two different SP's. The shape of the spectrum and the angle dependence of the gap have 
a striking similarity with the features observed in the normal state of the underdoped high-T c cuprates. 
We discuss also details about the phase diagram and the behaviour of the density of states. 

PACS. 7 4.25. -q 74.72.-h 74.25.Dw 74.25.Ha 



Many experiments performed for high T c cuprates pro- 
vide an evidence for the existence of a pseudogap in the 
underdoped regime above T c and below some temperature 
T* (5) which value increases with increasingthe doping dis- 
tance from the optimal doping, S opt —5 jp-[p|. The pseudo- 
gap is observed directly by angle-resolved photoemission 
spectroscopy (ARPES) measurements The strik- 

in g a bout this gap is its increase with increasing S opt — 
S [ |12fl while the critical temperature of superconducting 
(SC) transition, T sc , decreases. Another prominent fea- 
ture is the so-called (7r, 0) feature discovered by ARPES: 
the electron spectrum around the saddle-point (SP) is flat 
and disappears above some threshold value of wavevec- 
tor J9|. Several hypothesis exist about possible origin of 
the pseudogap In this paper we present another 

explanation of this phe nomenon in the framework of the 
model developed in p!^|-|po|. In these works the concept of 
the Electronic Topological Transition in 2D system is de- 
veloped and applied for the explanation of various effects 
experimentally observed in High - T c cuprates. 

In the present paper we consider various ordered states 
appearing in the vicinity of ETT point. We show that the 
ordered "excitonic" phase 19 formed in a proximity of 
QCP in a 2D fermion system on a square lattice is char- 
acterized (on one side of it, 5 < 8 C ) by the electron spec- 
trum strikingly similar to that observed in the underdoped 
cuprates. The mentioned QCP corresponds to the electron 
concentration n c — 1 — 5 C at which Fermi level (FL) crosses 



saddle points (SP) in the bare spectrum. As shown in p9| , 
in the case of hopping between more than nearest neigh- 
bors this point is a point of a fundamental electronic topo- 
logical transition (ETT) for which singularities in thermo- 
dynamical properties (and the logarithmic divergence in 
density of states at u = 0) is only first quite trivial as- 
pect (related to the local change in topology of the Fermi 
surface (FS)). The other aspects of criticality related to 
the mutual change in topology of FS near two different 
SP's lead to a very asymmetric behaviour of the noninter- 
acting and interacting system on two sides of ETT being 
quite anomalous on the side S < S c . On the other hand, 
for realistic for the high-T c cuprates ratios of hopping pa- 
rameters t' /t, 6 C is given by : S c — 0.27 for t'/t — —0.3 and 
S c = 0.17 for t'/t = —0.2, i.e. the anomalous regime S < 5 C 
occurs in the doping range where the experimentally ob- 
served "strange metal" behaviour takes place. Moreover, 
S = S c corresponds to a maximum of T sc (5) (as discussed 
in ]lSt| ) and therefore the latter regime can be considered 
as an underdoped regime. 

Some anomalies concerning the ordered "excitonic" 
phase have been discussed in 18 - j2^] . Namely, it was 
shown that the line of the "excitonic" instability grows 
from QCP to the side 6 < S c instead of having the form of 
a bell around QCP as it usually happens for an ordinary 
QCP. Other anomalies which exists in the ordered phase 
are considered in the present paper. [We call this phase 
"excitonic" ordered phase because the discussed instabil- 
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ity has the same origin as the classical "excitonic" insta- 
bility intensively discussed in the 60-70 p2|-]27|. Namely 
it is related to the opposite curvature of two parts of elec- 
tron spectrum in a proximity of FL. In the case considered 
they correspond to spectra in vicinities of two SP's.] 

We consider various possibilities for the ordered state, 
namely, Spin and Charge Density Wave orderings with dif- 
ferent types of the order parameter symmetries (s-wave, 
d-wave) depending on the effective interaction between 
the quasiparticles. Despite of the different symmetries, the 
properties of such ordered states resemble an "excitonic" 
states [p2[-[p7| are quite similar. The explanation of this 
phenomenon is in the fundamental role of ETT and ef- 
fects of criticality in the vicinity of the corresponding QCP 
poll . We show that the electron spectrum in the ordered 
phase is characterized by a gap on FL which opens at 
"hot spot" and extends until SP whatever is the doping 
concentration. Therefore there is always a gap at the SP 
wavevectors (0,±7r), (±7r,0). This remarkable feature is 
related, as we show in the paper, to a quite nontrivial as- 
pect of ETT: it is the end point of two critical lines for the 
"polarization operator" characterizing a behaviour of the 
free electron system. The other side of the same effect is 
an increase of the amplitude of the order parameter (and 
of the gap) with increasing the distance from QCP on the 
underdoped side. We show also that the electron spec- 
trum in a vicinity of SP gets a specific " flat" form which 
on one hand is typical for an " excitonic" phase (see for ex- 
ample ]2q| ) being a result of a hybridization of two parts 
of the bare spectrum with the opposite curvature and on 
the other hand has a striking similarit y w ith the form of 
the spectrum observed by ARPES We show that 

the spectrum " disappears" above some threshold value of 
wavevector in the direction (tt, 0) — (tt,tt) that is also an 
effect of the same hybridization. We briefly discuss also 
features related to strong-coupling limit of the model and 
effects of strong electron correlations. 

A starting point is a 2D system of free fermions on a 
square lattice with hopping between nearest (t) and next 
nearest (f) neighbors 



ek = — 2t(cos k x + cos k y ) — At' cos k x cos k y (1) 

(as in (T^l we consider t > and t' < 0) and \t' /t\ not 
too small and not too close to the limit = .5. The 

dispersion law ([[]) is characterized by two different saddle 
points (SP's) located at (± tt, 0) and (0, ±7r) (in the first 
Brillouin zone (— 7r,0) is equivalent to (tt, 0) and (0, — tt) 
is equivalent to (0, 7r)) with the energy e s = At' . When we 
vary the chemical potential fi or the energy distance from 
the SP, Z, determined as 



Z = fi- 



At' 



(2) 



the topology of the Fermi surface changes when Z goes 
from Z < to Z > through the critical value Z = 0, see 

In |19fl we have shown that such a system undergoes 
a fundamental ETT at the electron concentration corre- 
sponding to Z = Q. The corresponding quantum critical 



point is quite rich. It combines several aspects of criti- 
cality. The first standard one is related to singularities in 
thermodynamic properties, in density of states at u) = 
(Van Hove singularity), to additional singularity in the 
superconducting (SC) response function, all reflect a local 
change in the topology of FS. This aspect is not impor- 
tant for the properties we are interested in the present 
paper. Important aspects which reflect a mutual change 
in the topology of FS in the vicinities of two SP's are the 
following. First of all, it is a logarithmic divergence of the 
polarizalibility of noninteracting electrons 



°n, , A - 1 V" nF ( g q) ~ nF ( g q+k) 



as k = Q = (tt, tt) , lu = and Z — > 



X°(Q,0)cxln' 



(3) 



(4) 



which has an "excitonic" origin (uj max ~ t is a cutoff en- 
ergy). By "excitonic" origin we mean that two branches 
of the spectrum corresponding to vicinities of two SP's 
(a = t-2t', b = t + 2t') 

ei(k) = ei(k) — [i = —Z + ak x — bky, 

i 2 (k) = e 2 (k)-H=-Z + ak y '-bk 2 x (5) 

have such a form (see Fig.||) that at Z = the chemical 
potential lies on the bottom of one "band" and on the top 
of the another for the given directions (0, tt) — (tt, tt) and 
(tt, 0) — (0, 0), (see Fig.pl). Therefore, no energy is needed to 
excite the electron- hole pair. It is this divergence that is at 
the origin of density wave (DW) instability. The DW insta- 
bility can be of Spin Density Wave (SDW), Charge Den- 
sity Wave (CDW), Spin Current Density Wave (SCDWjor 
Orbital Current Density Wave (OCDW) instability Q) 
of interacting electron system depending on a nature of 
interaction. 

The nontriviality stems from the aspect of criticality 
related to the effect of Kohn singularity in 2D system : the 
point Z = 0, T = is the end of the critical line Z < 
each point of which is a point of static Kohn singularities 
in polarizability of noninteracting electrons. As shown in 

j, the latter aspect is a motor for the anomalous be- 
haviour of the system on the side Z > of ETT. One 
among the anomalies found in [|l9) concerns the ordered 
DW phases. We have obtained that the line of DW "ex- 
citonic" instability Tjjw(Z) has the anomalous form on 
the side Z > : it grows from QCP instead of having the 
form of a bell around QCP as it usually happens in the 
case of ordinary QCP and as it indeed happens on the side 
Z < 0. Below we show that this latter aspect is also at 
the origin of anomalous behaviour of the order parameter 
and of some other anomalies in the ordered state in the 
same regime Z > 0. 

on the side Z > of the electronic 



As shown in |19 



topological transition, a maximum of the static electron- 
hole susceptibility occurs at the wavevector q = Q . There- 
fore in a presence of q independent interaction or q depen- 
dent interaction negative for q ~ Q, the DW instability 
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happens at q = Q and this is the wavevector of ordering 
in the DW phase. As usual for such phases, one should 
consider a matrix electron Green function containing as 
components the normal and anomalous Green functions 
in terms of operators a£ and a\^ a which are the creation 
and annihilation electron's operators respectively: 



KuQL,iu n ) = - / dre^ T < T T a k(T (r)| a+ a (0) > 
Jo 

K 22 {k,iu> n ) = - / dre^ T < r T a k+Q(T (r)|a+ (0) > 
Jo 



Kg (k,tw„) 



dre iu n r < T T a k+QCT (r)|a+ (T ,(0) > 



. ( 6 ) 

[Below we will omit spin indices in the Green functions 
keeping in mind that K 12 = K^ a for CDW and OCDW 
states and K 12 = K$ for SDW and SCDW states.] 

If the anomalous Green function Ki 2 is nonzero (that 
should be found selfconsistently) the explicit expressions 
for the two Green functions are as follows 



u 2 (k) v 2 {k) 



lLO n - E\ UO n - E 2 

K 22 {k,ioj n )= + - 

_lUJ n — e 2 lUJ n - El 

K 12 (k,iu n ) = K 21 (k,iuj n ) = 

= u(k)«(k) 

[iuj n ~ei iuj n ~e 2 

where u, v - coefficients have a standard form 



(7) 



« 2 (k) = r 



! (k) = 



e A {k) - e B (k) 
2E{k) 

_ e A {k) - e B {k) 
2E{k) 




E{k) 



The spectrum in the ordered state is given by 



(8) 



CA + £b . ( e A -e B 

£1,2 = ^ ± \ ^ 



\A(kW 



e A (k) = e(k) e B (k) =e(k + Q), (9) 
where e(k) is defined by (Q). The equation for the gap is 

A(k) = -Tj2j^J2 ^(k' k + Q: P)^12(P, ^n) (10) 



where 7"i2 is a vertex which in mean field approximation 
coincides with the bare interaction: 

A 2 (k,k + Q,p) = Vq, forSDW{CDW) 
A 2 (k,k + Q,p) = Vk_ p , for OCDW {SCDW) (11) 

where Vk = 2V (cos(k x ) + cos(k y )). The type of the inter- 
action and therefore, type of the excitonic phase depend 
on the model. The SDW and OCDW instabilities occur 
in the case of a positive interaction in the triplet channel 
(exchange interaction), the CDW and SCDW instabilities 
take place for positive interaction in the singlet channel 
(density-density interaction). We will not fix for the mo- 
ment a type of interaction and therefore a nature of the 
ordered phase assuming that there exists either the first 
or the second interaction. 

The equation ( |l0| ) is reduced to the following equation 

l = 4\V\n°J (Q,Z,A) (12) 

where the " polarization operators" II DW (Q, Z, A) are given 
by one of the following equations [BO : 



n 



SDWfiDW 



(Q,Z,A) 



E-— — [tanh( 



47V 



2T 1 



tanh( 



-) 

2T' 



(13) 



77 



OCDW, SCDW 



(Q,z,A) = 



1 



(cos Px — cos Py ) e± 
47V ^ 475 p 1 y 2T' 



tanh( 



-) 

2T' 



_ (14) 
The expressions (fi~2) are the equation for the SDW, CDW, 
OCDW or SCDW gap which should be solved selfconsis- 
tently. We emphasis that for V > only SDW (OCDW) 
solution is possible whereas CDW (SCDW) solution takes 
place for V < 0. 

The solution of (|T2])-(14) is given by one of the follow- 
ing expressions 



A = A 



SDW 



(k) = A* 



SDW 



A = A CDW (k) = A$ DW , 

A = A OCDW {k) = A° CDW (cosk x - cosfc^/2, 

A = A SCDW (k) = A$ CDW (co8k x - cosfc y )/2 (15) 

The equations (^0|)-(|l5|) are quite standard. A non- 
triviality, as we show below, is related to the behaviour 
of the "polarization operator" in a proximity of ETT. 
As we have shown in |t9|, the effect that the point of 
ETT is the end point of the critical line Z < leads to 
the anomalous behaviour of the electron-hole susceptibil- 
ity x°{Q,Z 1 uj) on the side Z > 0. Below we show that 
a similar effect takes place for the "polarization opera- 
tor" (|l4|). The two functions coincides in the limit cases: 
X {Q,Z,uj = 0) = n DW {Q,Z,A = 0). [It is important to 
emphasize that the behaviour of the " polarization opera- 
tor" depends only on properties of the system of nonin- 
teracting electrons, namely on the topology of FS.[ 
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Calculated for T = "polarization operators" 7J DW/ 
(Q, Z, Aq) as a function of Ao for fixed Z (in the regime 
Z > 0) are shown in Figj|[ Since the properties of the "po- 
larization operators" are similar in many aspects we shall 
omit later the indices (SDW, CDW, OCDW or SCDW) 
except for the cases when it will be necessary to empha- 
size the difference. 

One can see that there is a singularity at some point 
Ao = A C (Z). The value of A C (Z) increases with increasing 
Z. The situation is quite similar to that analyzed in p9| 
for x° as a function of oj for fixed Z and T — 0. In the 
latter case we have found a square-root singularity at 

2Z 



1 - 2t'jt' 



which is the dynamic Kohn singularity. As we see in Fig. 3, 
for the polarization operator 7T(Q, Z, A) the singularity is 
weaker, while A C (Z) also scales with Z . 

Analytical estimations show that A C (Z) is given by 



A C (Z) = Z 



(16) 



while the asymptotic form of n(Q, Z, A) near the singu- 
larity is given by : 



tn{Q, Z, A) = 



_ / A X \\-AIA C 



A 2 \l-A/A c 



B, 
B, 



A<A C 

A > A c . [U) 



The jump in the derivative, A\ — A 2 , depends only on t' jt 
]3l| and is proportional to 



Ax - A 2 oc 



1 



In I 



A 



(18) 



1*7*1 V *'/* 

where Ao is a constant (for the spectrum (||) Ao ~ 7r/8). 
The critical line (|l6|) is clearly seen in Fig. 4a where we 
present the calculated 27(Q, Z, Aq) as a function of Z and 
Ao- From the point of view of the behaviour of the "po- 
larization operator" , the ETT point is the end point of 
two critical lines. The first is the semiaxis Z < each 
point of which corresponds to the square-root singularity 
in n(q, Z, A) occurring as A — ► and q — > q m , where 
the latter is the characteristic for this regime wavevector 
of incommensurability (see 1 19 where the q dependence of 



-Z7(q, Z, 0)= x (q 5 Z) is analyzed in details.) The second 
is the line Z — Ao each point of which corresponds to the 
kink in 77(q, Z, A(q)) occurring at T — as q — » Q where 
the latter is the characteristic wavevector for the regime 
Z > 0. At the point of intersection of these lines, Z = 0, 
the two types of singularities are transforming into the log- 
arithmic singularity ; 7T(q, 0, A) oc In | max(q — Q, A)\. 

The existence of the growing with Z critical line de- 
termines a quite unusual form of the lines II(Q, Z, A) = 
constant which develop around the critical line A C {Z) and 
grow with increasing Z, (see Fig.||(b)). 

In preceeding discussion we presented some general 
analysis which does not depend on details of interaction 
considered but only on the topology of the FS. To pro- 
vide the calculations, let us consider a particular case of 



interaction resulting in Spin Density Wave (pd|) , (|l3j) . The 
solution of corresponding eq. ( |l2| ) for t/V — 1.8 is shown 
in Fig|| Two branches of the solution have an anomalous 
dependence of the gap on Z reproducing the form of the 
lines II(Q, Z, Ao) = constant in Fig.[|(b). The anomaly is 
that for both solutions gap increases with increasing 
the distance from the quantum critical point, i.e. 
from the point which is at the origin of the ordered phase. 
[For an ordinary QCP the gap is maximum at the elec- 
tron concentration corresponding to QCP and decreases 
monotonously with increasing the distance from QCP. For 
example such a picture takes place for DW phase on both 
sides from QCP in the case oft 1 =t" = ... = 0; as we dis- 
cussed in ]1!| in the latter case all anomalies in the regime 
5 < 8 C disappear. In the case considered in the paper it 
happens on the overdoped side of the QCP.] 

The difference between two solutions for the gap pre- 
sented in FigJ| is that 



while 



Z\i(Z) > z 



A 2 (Z) < Z 



(19) 



(20) 



for any Z, any t/V , any t' jt since the two lines, A\(Z) and 
A 2 (Z) are attached to the critical line A = A C (Z) = Z 
from above and from below. For the most range of the ex- 
istence of the ordered phase Z < Z$, see Fig. 5l only one 



solution exists, the one corresponding to eq. 



In the 



hyperbolic approximation and under the condition \t' /t\ 
not too small Z^ is given by: 

Zg> cxw ma:c exp(-7r 2 t/(yin|<A , |)). 



For this solution one has 

A 1 (Z) = Ao(Z)=f(Z) + A(0) 
where A(0) is given by 



A(0) oc |i'|exp(-- 



2TT 2 \t'/V\ 



f), 



(21) 



(22) 



and f(Z) is an increasing function of Z, linear under the 
condition A(Z) 3> A(0). The expression (|22| ) is valid un- 
der condition ir 2 \t' /V\/ -y/l - (2t'/t) 2 > 1. For the narrow 
Z range of the coexistence of the two solutions Z^ < 
Z < Zc 2 ) it is the solution Ai which is favorable (see Ap- 
pendix). Therefore, the value of the gap increases 
with increasing Z being always larger than Z. As 
we have shown, this is a consequence of the effect that the 
point of ETT is the end point of two critical lines. 

Let's analyze now the form of the spectrum in the DW 
phase. The spectrum given by (||) is plotted in Fig.^. for 
three important directions : (n, ir) — (ir,0) — (0,0) and 
(0, 0) — (it, it). The spectrum in the vicinity of SP has the 
following prominent features: The first is a characteristic 
" flat" shape (very close to the experimental shape Q , see 
Fig.||(a) being a consequence of the hybridization of the 
two branches of the bare spectrum in the vicinity of two 
different SP's with the opposite curvatures, (see Fig.|^). 
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The second: the spectrum in the direction (it, tt) — (0, tt) 
"disappears" above some threshold value of wavevector 
since the residue u£ tends to zero (that is also an or- 
dinary consequence of the hybridization). On the other 
hand, since 

ei(ksp) = -Z + A, 

e 2 (k S p) = -Z-A 

(see, e.g. Eq.(||)) and A > Z, the chemical potential 
always lies in the gap for the part of Brillouin zone 
(BZ) starting from the "hot spot" until SP that is a 
consequence of the existence of the critical line A — A c 
related to the discussed above aspect of criticality of the 
QCP. For the direction (0,0) — (1, 1), Fermi level lies in 
the lower branch of the spectrum, (see Fig.||(b)), i.e. the 
system remains metallic. The theoretical spectrum has a 
striking similarity with the anomalous experimental elec- 
tron spectrum in the underdoped cuprates observed below 
the characteristic line T* (5) by ARPES || , we reproduce 
it in Fig.||. We remind that ARPES measures a spectral 
function only below FL. 

Then in Fig.^ we present the angle dependence of the 
value of £k — i-e. of the gap calculated from FL, in 
the same way as it is done in ARPES experiments ]l0| . 
Namely we plot the minimal value of |ek — /j, for each 
given direction. The dependence is of a "d-wave type" in 
a sense that the gap increases with increasing the argu- 
ment (cosfc x — cosfcj,) almost linearly in the proximity of 
SP. However the dependence is flat (not linear as it hap- 
pens in the d-wave case) when approaching the direction 
(1,1). Such a behaviour is also close to the experimentally 
found behaviour above T c Il0| reproduced in Fig.||(b). [Al- 
though the authors of |10|| claim that the behaviour ob- 
served above and below T c is the same, what one sees in 
the experimental plot is not exactly this : the behaviour 
above and below T c is similar in the vicinity of SP and 
different when approaching the (1,1) direction and this 
occurs quite systematically, see also the plots in jl0| for 
other samples.] 

We considered the particular case of SDW as an exam- 
ple of ordreed "excitonic" state. Nevertheless, all aforesaid 
is true for any other types of ordered states since the ex- 
istence of such states is determined only by topology of 
FS. 

Let us study now the one particle density of states 
(DOS) given by the expression 

^ = --^E[^n(P> £ ) +IrnK*(p,e)] = 



N 



£ - £i(p)) + S(e - £ 2 (p))], 



(23) 



Numerical calculations with the spectrum (|l]) give the 
picture shown in Fig.[To|. 

The density of states of SDW (CDW) states deviates 
from the DOS in the initial metallic state in two e ranges 
notated as A and B. For OCDW (SCDW) states only 
feature A survives. Analytical calculations show that the 



A-feature is related to the existence of the discussed above 
QCP (which we call below QCP1). Calculations of the 
integral in ( p3| ) performing with the hyperbolic spectrum 
(||) valid in the vicinities of SP's show that in the A range 
DOS is characterized by three singularities (instead of one 
logarithmic singularity in the bare density of states ATq(s) 
as e — > — Z). Those are a logarithmic singularity at E\ = 

-z-A Q ^i-A(t>/ty m 



M(e->ex) 



dn(- 



(24) 



and jumps at two energies 

£ 2 ,3 = -Z±A Q . 

The distance between two jumps is equal to 2A. 

The B feature is related to the existence of the second 
quantum critical point in the system (QCP2) discussed in 
[ p9[ . This point corresponds to the electron concentration 
when the chemical potential is equal : ^ = fi C 2 = or 
by other words when the wavevector connecting two parts 
of FS in the direction (1,1) is equal to Qaf — {tt,tt). 
In this case two "hot spots" on FS come together at the 
singular position (±7r/2, ±7r/2) before disappearing. The 
calculations of the integral in ( J23| ) with the spectrum taken 
around (7r/2,7r/2) give a logarithmic divergence at the 
point £ 4 = -Z - At' + A(ir/2, tt/2): 



N{e £ 4 )-A/- (£4) ~ -J ^ m( F - 7Ij 



) (25) 



and a jump at the point e = —Z — At' — A(n/2, tt/2). This 
feature does not exist for OCDW (SCDW) states since 
A(k) = along the diagonal of BZ. 

The B feature is important in the case when the chem- 
ical potential lies close to the pseudogap in the B part that 
should take place in the electron-doped cuprates. For the 
hole-doped cuprates we are interested in the present pa- 
per, it is QCP1 which determines properties of the system. 
In this case the chemical potential lies in the "pseudo-gap" 
A according to the properties of the electron spectrum in 
the vicinity of SP discussed above. 

Let's analyze now the range of the existence of the or- 
dered phase in the T — Z plane. For this sake let's analyze 
the behaviour of 7T(Q, Z, Aq) as a function of Z at finite 
temperature. [We again consider SDW state for certain- 
ity.] Results of calculations are presented in Fig|ll]. The 
first observation is that the gap changes only little with 
T at low T. The second is that the behaviour at finite 
temperature as a function of Z is qualitatively the same 
as for T ~ and it is anomalous: the value of the gap 
increases with increasing Z. 

The phase diagram in T — Z plane obtained for SDW 
(CDW) instability based on the analysis of the gap be- 
haviour at finite T is presented in Fig. 12. It is worthwhile 
to note that the polarization operator 77 (Q, Z, A) ( p^ ) 
calculated for OCDW (SCDW) ordered states has essen- 
tially more abrupt behavior as a function of Z in com- 
parison with those for n SDW - CDW {Q, ZA) (01). Such be- 



haiour appears due to additional factor (cos(p x )— cos{p y )y 
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in the integral (U3). As a result, the domain of existence 



of OCDW (SCDW) solutions for equation fll2|) at various 
doping concentrations is substantially narrower than for 
SDW (CDW) case. Nevertheless, it does not affect on the 
qualitative shape of phase diagram Fig.|l2"[ 

The solid line is the line where A\{T) = 0. The dashed 
line is the line where A$(T) = 0. These two lines are at 
the same time the lines of instabilities of the undistorted 
metallic state. The line A^{T) = is not however a line 
of a phase transition since the nonzero solutions for the 
gap exist on the left of this line until the dot-dashed line. 
Along the latter line corresponding to the disappearance 
of the "ordered" solution, the gap is finite and the two 
solutions coincide : A (T) = A^T) = A 2 (T). The situa- 
tion is clearly seen from Fig jl^ where we present the lines 
7T(Q, Z, A) = const for different Z and fixed tj J which in 
fact give the full picture of the behaviour of the DW gap 
as a function of Z and T. 

As we discuss in the Appendix, in the region between 
the dot-dashed and dashed line, where three solutions 
Aq = A\, Aq = Ai and A — coexist, it is the solu- 
tion Z\o = A\ which is energetically favorable. 

Thus, the dot-dashed line in the phase diagram in 
Fig.[l^ is the line of the first-order phase transition. The 
gap along this line changes only little at low temperature 
and tends to zero rapidly in the vicinity of the point O. 
The latter is a tricritical point. The range in T — Z plane 
in the vicinity of this point corresponds to a strongly fluc- 
tuating regime which we will consider elsewhere. It is im- 

(2) 

portant to add also that at the point Z = Z cr ' of the 
appearance of the ordered phase at T = 0, the gap is ex- 
actly equal to Z that means that the upper branch of the 
spectrum in Fig.|6|(a) touches FL. Then when moving in- 
side the ordered phase the gap A becomes larger than Z 
and this branch goes up leaving the FL. 

Above we have considered the critical temperatures 
and the gap behaviour as functions of the energy distance 
from the QCP, Z. It is worth for applications to cuprates 
to change the description and to consider physical prop- 
erties as functions of electron concentration n e or of hole 
doping 5 = 1 — n e . To do this we use the relation between 
Z (or the chemical potential a) and the hole doping which 
for T — is given by : 



1 



N(uj)duj. 



So far as 



Z oc 5 C 



(26) 



(27) 



all dependencies considered above can be rewritten as 
functions of doping distance from QCP. For example, the 
phase diagram in the plane T— 6 calculated for t' ' Jt = —0.3 
for which S c = 0.27 gets the form shown in Fig .[14]. 

One can easily obtain values of doping for all plots pre- 
sented in Figj4j-Fig.|ll"| when comparing the phase diagram 
in T — Z plane in Figjl^, and in T — 5 plane in Fig.|l4|. 

Obviously, the gap Aq(S) increases with 8 C — 5 in the 
same way as it increases with Z, see Fig|^ and Figju] for 
A = A x . 



All features discussed above do not depend on the na- 
ture of the ordered phase, SDW, CDW, OCDW or SCDW 
since they reflect the topological aspects of ETT. The type 
of the excitopic phases developing around ETT point de- 
pend on the type of interaction. It is the SDW or OCDW 
state in the case of a positive interaction in the triplet 
channel (exchange interaction) and the CDW or SCDW 
state in the case of a positive interaction in a singlet 
channel (density-density interaction). The ordered SDW 
phase is characterized by spin ordering with momentum 
< S 7 q >= l/2(< n aa {Q) > - < n ss (Q) >) = A and 
the CDW phase by the charge ordering. In the SCDW 
(OCDW) the staggered magnetization (density) is equal 
to zero. Nevertheless, the spin-current (charge-current) 
correlation functions survive. 

In our opinion for the case of high-T c cuprates it is the 
interaction in the triplet channel which determines the be- 
haviour of the system and the nature of DW phase. From 
the theoretical point of view it is this situation which cor- 
responds to the strong-coupling limit models : the Hub- 
bard model and the t — J model. For example for the 
latter with the J term written as Hj = Y)^ Jij{aSiSj — 

(b/^rnrij} one has V^ DW = aJ q while V° DW = -|j q , 
i.e. the interaction in the triplet channel is positive while in 
the singlet channel is negative. This version is supported 
also by experiments in the high-T c cuprates : observed 
experimentally (by INS, see for example |33| and NMR) 
strong magnetic response around q = Q is a phenomeno- 
logical argument in a favor of a strong momentum de- 
pendent interaction in a triplet channel, i.e. of Vq = J q 
(J > 0). However, we can not exclude an importance of 
an interaction leading to the CDW (SCDW) order. 

Another point concerning the interaction is its strength. 
Depending on the ratio |V|/W (where W is an energy 
bandwidth), maximal T^$j can be high or low. Respec- 
tively, the DW phase can lean out of SC state or can be 
hidden under it. [In the presence of the interaction in the 
triplet channel, J q , both SDW and SC instabilities occur 
around QCP1 under the same condition : J > 0, for the SC 
instability see It is temptating to identify the prop- 
erties obtained for the DW state with the properties ob- 
served experimentally in the underdoped cuprates above 
T sc (6) and below T*(S). Indeed they have a striking resem- 
blance, as one can see when comparing FigJ^ and FigJ^, 
Fig.^(a) and Fig.^(b) and when comparing the behaviour 
of the gap as a function of Z (or doping, S c — S) with 
the experimental behaviour JlJ]. Our calculations (when 
considering both d-wave SC and DW instabilities in the 
presence of interaction J in the triplet channel) show that 
the answer is quite subtle. When t' jt = —0.2 the ordered 
DW phase leans out of the SC phase for tj J < 1.90, for 
t'/t = —0.3 this happens when t/J< 1.55. So far as re- 
alistic value of tj J for cuprates is estimated to be in the 
interval t/J = 1 — 3, both variants (when the DW phase 
leans out or is hidden under the SC phase) are possible. 

In this case two scenarios can be discussed. First, the 
DW state of s- or d- wave symmetry can coexist with d- 
wave SC state (this situation is considered in |36|) result- 
ing in appearance of supplementary 7r-triplet state. Sec- 
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ond, the DW state can be suppressed under the SC state. 
Nevertheless, the strong "fluctuation" memory of this or- 
dering will affect on the behaviour of disordered metallic 
state. 

Even in the case if the long-range ordered DW phase is 
hidden under SC phase it is this type of ordering which de- 
termines short range correlations in the disordered metal- 
lic state above T c and below T*(S). This point is discussed 
in Iplhlpj. The state below T*(S) (x5 c -S and above T c is 
quite exotic. It is almost frozen in both temperature and 
doping. By this we mean that the parameter k 2 which 
determines a proximity to the ordered DW phase does al- 
most not change neither with T nor with doping remaining 
therefore quite low in a wide region in T — 5 plane below 
T*(S). Such a quasiordered state keeps a strong memory 
about the ordered phase. Therefore, electron properties in 
this state should be close to those in the DW ordered state 
being however characterized by strong damping. [By the 
way it is exactly what is observed by ARPES. The ex- 
perimental electron spectrum has a form shown in Fig.^|. 
being however characterized by a spectral function of a 
very damped form. Explicit consideration of the electron 
spectrum in this state will be presented elsewhere]. 

Summarizing, we have studied the DW phase which is 
formed around QCP1 (associated with ETT) and we have 
shown that this phase is characterized by the following 
prominent features: 

(i) the specific "flat" shape of the spectrum in the 
vicinity of SP, 

(ii) " disappearance" of the spectrum above some thresh- 
old value of wavevector in the direction (jr, 0) - (jr, 7r), 

(iii) pseudogap in DOS with FL lying in it, 

(iv) increasing of the gap in the spectrum around SP 
and of the pseudogap in DOS with decreasing doping for 
5 c -5 

(v) angle dependence of the gap calculating from FL 
which is of a d-wave type close to SP and flat close to the 
direction (1,1). 

All these features have a striking similarity with the 
experimental features revealed by ARPES in the normal 
state of the underdoped hole-doped cuprates. 



Appendix A 

The free energy density in the approximation correspond- 
ing to considered in the paper is given by: 



F =~ ^EEH2cosh(£^l)) + f] 



k a=l,2 



fj,n. 



(28) 



[Note that the equation ( |i"o|) corresponds to dF/dA = 
Therefore, the difference between free energies correspond 
ing to A = A\ and A = Ai is given by 



F-i —Fi 



A\-A\ 
W 



1 cosh(- 



k a=l,2 



2T 



^ cosh( 



g a (k, A 2 

2T 



One can check by numerical calculations that F± — 
F 2 < for the whole range of the coexistence of the two 
solutions. 

Some analytical estimations can be also done for low T 
based on the well-known expression |35| for the difference 
between thermodynamic potentials of the ordered and dis- 
ordered states: 



5Q = Q{Ax) - J7(0) = / 
Jo 



h d(l/V) 
dA 



Ad A. 



(30) 



When substituting the expressions for A\ (pl|), (22) one 
gets 



SF/t = 6f2/t 



v / i-4|t'AI 2 4? 



\t'/t\ t 2 A{Z = Q) t 2 A(Z = 0)' 

(31) ' 

One can see that this correction is negative. Therefore, the 
solution A = Ax is favorable with respect to the solution 
A = for any < Z < z£ ] . 

* Present address: Institiit fur Theoretische Physik, Uni- 
versitdt Wiirzburg, D-97074 Wiirzburg, Germany 
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Fig. 1. Fermi surface of the electron system with the dispersion law ([!]) for different Z and t' /t — —0.3. The thick line 
corresponds to Z — 0. 




Fig. 2. Schematic presentation of the electron spectrum in a vicinity of two SP's for Z — 
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Fig. 3. Calculated "polarization operator" 77(Q, Z, A(Q)) as a function of zAq for fixed Z and T = 0. 
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Fig. 5. Gap A obtained by solving equations (|l2|), ( |14| ) as function of Z (t/V=1.8, t' /t = —0.3). The solid line corresponds 
A\(Z), the dot-dashed line to A2(Z), the dashed line to A C (Z). 
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a) b) 

Fig. 6. One particle spectra along (tt, tt) — (n, q y /n) and (tt — q x /^, 0) — (0, 0) symmetry lines (a) and in (0, 0) -(tt, tt) direction 
(b), t'/t = — 0.3, t/V = 1.8, Z/t — 0.03. Long dashed line is the bare spectrum, dot-dashed line corresponds to the spectrum 
when the residue of the Green function (^) less than 0.1. 




Fig. 7. Schematic representation of the bare spectrum in the vicinity of the two saddle points for 
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Fig. 9. Theoretical angle dependence of the SDW gap calculated from FL in the underdoped regime Z > (t' jt — —0.3, 
t/V = 1.7, Z/t = 0.3) (a) and the experimental leading edge midpoint measured by ARPES in the underdoped BSCO [11] (b). 
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Fig. 10. Density of states in the ordered "excitonic" phase calculated for Z/t = 0.03 (t/V = 1.8, t'/t = —0.3). Dashed line 
corresponds to the DOS in the initial metallic state. 
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Fig. 11. The DW gap in t units as a function of Z for increasing temperature : T/t = 0.005, 0.1, 0.2 (t'/t = -0.3, t/V = 1.8). 
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0.06 




Fig. 12. Phase diagram around QCP1 in T-Z coordinates (t/V = 1.8, t'/t = —0.3). We show only the regime Z > 
corresponding to the anomalous behavior. The solid line is a line of second-order phase transition, the dot-dashed is a line of 
first-order phase transition and the dashed line is a line of instability of the disordered metal state (spinodal). The point O is 
a tricritical point. 




Fig. 13. Lines of 7T(Q, Z, A) = const for fixed t/V = 1.8 and different Z. The plot (b) is a zoom of the plot (a) corresponding 
to the coexistence of the two solutions for the gap. 
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Fig. 14. Phase diagram around QCP1 in T-S coordinates (t/V = 1.8, t' jt — 



